Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN MfeK style file vl.4) 



Relativistic Distortions in the X— ray Spectrum of Cyg X— 1 



C. Done 1 , P.T. Zycki 1 ' 2 

1 Department of Physics, University of Durham, South Road, Durham, DH1 3LE 

2 Nicolaus Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw 



5 February 2008 



00 

o 

Q 



(N 
(N 

oo 

ON 

a: 

6 
to 

C3 



13 



ABSTRACT 

We present the first significant detection of relativistic smearing of the X-ray 
reflection spectrum from the putative accretion disk in the low/hard state of Cyg 
X— 1. The ionization state, and amount of relativistic smearing are simultaneously 
constrained by the X-ray spectra, and we conclude that the disk is not strongly ionized, 
does not generally extend down to the last stable orbit at 3 Schwarzschild radii and 
covers rather less than half the sky as seen from the X-ray source. These results are 
consistent with a geometry where the optically thick disk truncates at a few tens 
of Schwarzschild radii, with the inner region occupied by the X-ray hot, optically 
thin(ish) plasma. Such a geometry is also inferred from previous studies of the reflected 
spectrum in Galactic Black Hole transient sources, and from detailed considerations 
of the overall continuum spectral shape, suggesting that this is a robust feature of 
low/hard state accretion onto Galactic Black Holes. 

Key words: accretion, accretion discs - binaries: close - black hole physics - stars: 
individual (Cygnus X-l) - X-ray: general - X-ray: stars 



1 INTRODUCTION 

Some of the strongest evidence for the existence of black 
holes has come from recent ASCA (0.6-10 keV) X-ray ob- 
servations of the shape of the iron Ka line in Active Galactic 
Nuclei (AGN). AGN typically produce copious hard X-ray 
emission, and the iron line is formed from fluorescence as 
these X-rays illuminate the infalling material. The combi- 
nation of Doppler effects from the high orbital velocities and 
strong gravity in the vicinity of a black hole gives the line a 
characteristically skewed, broad profile (Fabian et al., 1989). 
This has been unambiguously identified from ASCA 0.6-10 
keV spectra of the AGN MCG-6-30-15, where the breadth 
of the line profile implies that the accretion disk extends 
down to at least 6R S , the last stable orbit in a Schwarzschild 
metric (Tanaka et al 1995; Iwasawa et al 1996). 

As well as producing the line, some fraction of the illu- 
minating hard X-rays are reflected from the accretion flow, 
producing a characteristic continuum spectrum. The am- 
plitude of the line and reflected continuum depend on the 
amount of material being illuminated by the hard X-rays, its 
inclination, elemental abundances and ionization state (e.g. 
Lightman & White 1988; George & Fabian 1991; Matt, Per- 
ola & Piro 1991; Ross & Fabian 1993; Zycki & Czerny 1994). 
The mean amount of reflection and line seen in Seyfert 1 
AGN is consistent with a power law X-ray spectrum illumi- 
nating an optically thick, (nearly) neutral disk, which sub- 
tends a solid angle of ~ 2ir (Pounds et al 1990; Nandra 



& Pounds 1994), although it is now recognised that indi- 
vidual objects show significant dispersion about this mean 
(Zdziarski, Lubihski & Smith 1999). 

This contrasts with the situation in the Galactic Black 
Hole Candidates (GBHC). These are also thought to be pow- 
ered by accretion via a disk onto a black hole, and, in their 
low/hard state, show spectra which are rather similar to 
those from AGN. However, the amount of reflection and iron 
line is much less than would be expected from an accretion 
disk which subtends a solid angle of 2n (e.g. Gierlinski et 
al 1997a), and the detected line is narrow, with no obvious 
broad component (Marshall et al 1993, Ebisawa et al 1996, 
hereafter E96). There are (at least) two possible explana- 
tions for this: firstly that it is an artifact of the difference 
in ionization state of the disk between GBHC and AGN, or 
secondly that there is a real difference in geometry between 
the supermassive accreting black holes and the stellar mass 
ones. 

Ionization differences seems at first sight to be an ex- 
tremely attractive option. For accretion at the same fraction 
of Eddington, the disk temperature should scale as M _1//4 . 
The GBHC inner disk is then expected to be a factor of 
~ 30 hotter than in AGN, and the higher expected ioniza- 
tion state from the thermal ion populations gives a reflected 
continuum and associated iron line that can be very differ- 
ent to that from a neutral disk (Lightman & White 1988; 
Ross & Fabian 1993; Matt Fabian & Ross 1993ab; Zycki et 
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al 1994; Matt Fabian & Ross 1996; Ross, Fabian & Brandt 
1996). As well as increasing the energy of the iron edge (and 
line), ionization also reduces the low energy opacity of the 
disk, and so enhances its reflectivity. However, the iron edge 
(whatever its energy) has a fairly constant cross-section, so 
the net result is to increase the relative strength of the iron 
edge feature in the reflected continuum. There is also the 
possibility of the line being suppressed through resonant ab- 
sorption over a (rather small) range in ionization states, so 
the net result can be a prominent edge without obvious line 
emission accompanying it. Such models were proposed as the 
solution to the weak line emission seen in GBHC compared 
to AGN (Ross & Fabian 1993; Matt et al., 1993a; Ross et al 
1996). 

However, detailed spectral fitting with photo-ionized 
reflected continuua shows that the ionization state of the 
reflector is in general too low for the resonant absorption 
process to be important. This argues against such models, 
although relativistic shifts could affect the derived ionization 
states (Matt et al 1993a, Ross et al., 1996). Another way to 
suppress both line and reflected continuum is from a radial 
dependance of ionization state. If the inner region of the 
disk were so highly ionized that they produce no atomic 
features then their contribution to the reflected continuum 
is unobservable in the 2-20 keV range. Here we fit data 
from Cyg X-l with a reflection model which includes both 
ionization (as a function of radius) and relativistic smearing 
in order to test whether the spectra are indeed affected by 
resonant absorption and whether radial ionization structure 
is important in masking the reflection signature. 

We significantly detect relativistic smearing of the re- 
processed features in Cyg X-l, and find that the disk is not 
highly ionized. This apparently rules out both featureless 
(extreme ionization) reflection and resonant absorption. The 
line is not suppressed relative to the reflected continuum, it 
is merely broadened by the relativistic effects which made 
it difficult to detect in ASCA data (E96), where the ob- 
served weak and narrow line emission probably comes from 
the companion star. The derived covering fraction of the 
relativistic reflector is substantially less than expected from 
a flat disk. The most plausible explanation for this is that 
the optically thin(ish) X-ray corona fills a central 'hole', 
so that less than half of the X-ray flux intercepts the disk. 
The amount of relativistic smearing generally is inconsistent 
with the disk extending down to the last stable orbit at 6R g . 
These results from the shape of the reprocessed spectrum are 
remarkably similar to the geometry derived from energetic 
arguments from the continuum spectral shape (Gierlinski et 
al. 1997a; Dove et al 1997; Poutanen, Krolik & Ryde 1997). 
Equally, they are very similar to the derived reflection ge- 
ometry from spectra of transient GBHC during their decline 
(Zycki, Done & Smith 1997; Zycki, Done & Smith 1998ab), 
pointing to a single, fairly stable geometry for the low/hard 
state GBHC. 



2 THE SPECTRAL MODEL 

The reflection model is described in detail in Zycki et al 
(1998b). Basically it consists of the angle dependant reflec- 
tion code of Magdziarz & Zdziarski (1995), with ion popu- 
lations calculated as in Done et al (1992), as implemented 



in the 'pexriv' model in XSPEC. This determines the ion 
populations by balancing photoelectric absorption against 
radiative recombination. The photo-electric absorption edge 
energies for iron were corrected from the rather approximate 
values given by Reilman & Manson (1978) to those of Kaas- 
tra & Mewe (1994). As noted by E96, these give a significant 
difference in the derived ionization state. Photo-ionization 
rates are strongly dependent on the parameter £ = L/(nr 2 ) 
ergs s~ x cm -1 (where L is the ionizing luminosity, here taken 
to be between 5eV and 20 keV, n is the density and r the 
distance of the material from the ionizing flux) and weakly 
on spectral shape (assumed here to be a power law), while 
recombination depends on temperature (assumed fixed at 
10 6 K). The self-consistent iron line is then calculated from 
the Monte-Carlo code of Zycki & Czerny (1994) and added 
to the continuum. We assume fixed 'solar' abundances from 
Morrison & McCammon (1983), except for iron, which is 
free to vary from its tabulated value of Fe/H = 3.3 x 10~ 5 . 

Ionization as a function of radius can easily be treated 
by adding together many radial rings with ionization f (r) oc 
r 13 , where the relative covering fraction for each ring can be 
calculated from the X-ray irradiation Fi rr (r) oc r -3 . 

The relativistic smearing is taken from the XSPEC 
'diskline' model, with slight modifications to take light 
bending into account (Fabian et al 1989). This is param- 
eterised by the inner and outer radius of the disk, 7?j n and 
i?out, its inclination and (again) the irradiation emissivity, 
F irr (r) oc r~ 3 . We fix R out = 10 4 i? g (where R g = GM/c 2 ) 
and fit for i?i n , convolving the relativistic transfer function 
with the total (line plus continuum) reprocessed spectrum. 
For the multiple ionization state disk the relativistic transfer 
function is calculated for each annulus, convolved with the 
reprocessed spectrum from that ring and then co-added to 
get the total spectrum. 

We assume a power law illuminating spectrum. This is 
an approximation to the true shape of a Comptonised spec- 
trum and it becomes a poor description at energies close to 
the seed photon energy and close to the electron tempera- 
ture. Typically in low/hard state of Cyg X-l the electron 
temperature is ~ 100 keV (e.g. Gierlinski et al 1997a), so 
this will not give significant distortions to the spectral range 
considered here, which is always < 20 keV. 

The situation with regards the soft photons is less 
clearcut. The strongest component of the soft X-ray contin- 
uum can be modelled as a blackbody or disk blackbody with 
a temperature of 0.1 — 0.2 keV (e.g. Balucinska & Hasinger 
1991; E96; di Salvo et al 1998). Again, this is far enough 
away from the energy range considered in this paper (> 2 
keV) for the distortions of the true Compton scattered spec- 
trum from a power law to be negligible. However, there is 
also a second soft component which extends to higher en- 
ergies (E96, di Salvo et al 1998). E96 modelled this in the 
ASCA data using a softer power law below a break energy 
around 3 keV, but such a broken power law is rather un- 
physical. Perhaps a more realistic model is one where soft 
photons from the accretion disk (at ~ 0.1 keV) are Comp- 
tonised by a rather cool population of electrons with kT ~ 1 
keV as well as by the hot electrons with kT ~ 150 keV. The 
disk photons at ~ 0.1 keV are energetically more important, 
so it seems probable that this component at ~ 1 keV does 
not contribute significantly to the seed photons. 

The physical origin of the second soft excess is entirely 
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Figure 1. The ratio of each EXOSAT GSPC spectra of Cyg 
X-l to a best fitting power law model. Both excess counts at 
the line energy and a deficit at the edge are seen, as in AGN. 
The dotted line shows the same ratio for the GINGA-12 AGN 
spectrum (Pounds et al 1990) for comparison. 

unknown, and so is its geometry with respect to the disk and 
hot plasma. It may or may not then produce its own reflected 
component. We assume it does not, but note that even if it 
does then its contribution to the crucial iron line and edge 
region will be very small due to the low temperature of this 
component. 

In the following fitting we generally neglect data below 
4 keV so as to avoid the complexities of the soft excess, and 
fix the galactic column at 6 x 10 21 cm~ 2 . 



3 CONSTANT IONIZATION REFLECTION 
3.1 EXOSAT GSPC data 

The EXOSAT GSPC data from Cyg X-l give some of the 
best spectra to date from this object, with the broad 2-20 
keV energy range of GINGA data but resolution comparable 
to that of the ASCA GIS. We use the 5 GSPC spectra of 
Done et al., (1992). The residuals to the best fit simple power 
law for each GSPC spectrum are shown in Figure 1, together 
with the residuals to a simple power law fit to the GINGA- 
12 AGN spectrum of Pounds et al (1990). Contrary to claims 
by Ross et al (1996), Cyg X-l is not dominated by the edge 
structure, but also shows a strong excess at the iron line 
energy. Also, both line and edge appear equally diminished 
in Cyg X-l compared to AGN, although it is somewhat 
misleading to compare these residuals directly since Cyg X- 
1 is strongly absorbed by the interstellar medium and can 
have a substantial contribution at 2 keV from the observed 
soft excess. 

We fit the 5 GSPC spectra simultaneously in XSPEC, 
although only the reflection parameters of inclination and 
iron abundance are constrained to be equal across all the 
datasets. 

We first fit the 4-20 keV data with a power law. This 



gives a poor fit (xl = 1309/625) due to the presence of 
spectral features from iron (see Figure 1). Including a neu- 
tral reprocessed component inclined at 30° (Gies & Bolton 
1986) gives a dramatic decrease in xl to 677/620, with a 
further significant improvement in the fit if this is allowed 
to be ionized (xl = 611/615). Allowing the iron abundance 
to be a free parameter makes no significant difference, with 
xl = 610/614 for A Fc = 0.9 ± 0.2x solar. 

Allowing the reprocessed spectrum to be relativistically 
smeared gives a significantly better fit, with xl = 535/609. 
We obtain constraints on the inner radius of the disk in 
each spectrum, and these are somewhat larger than the 6R g 
expected for a disk extending down to the last stable or- 
bit, assuming an illumination which is cc r~ 3 . The derived 
iron abundance is Af c = 1.8j^J'.Ix solar i.e. supersolar abun- 
dances are preferred. This is consistent with the ~ 2x solar 
abundance of iron (given as Fe/H= 6 ± 1 x 10~ 5 ) inferred 
for the stellar wind of the companion star from the strong 
edge feature seen in the absorbed 'dip' spectra (Kitamoto 
et al., 1984). The relativistic smearing allows more iron line 
to be present, since it is broad rather than narrow, so it is 
much less observable. 

One potential problem may arise if there is a substan- 
tial narrow line component from X-ray illumination of flared 
outer regions of the disk and/or the companion star (Basko 
1978) and stellar wind. This could distort the relativistic line 
fit, by adding a sharp core to the line profile. Such a line is 
indeed detected in high resolution BBXRT and ASCA data 
(Marshall et al 1993; E96). We allow for this by including a 
narrow Gaussian line at 6.4 keV, whose normalisation is a 
free parameter for each spectrum. This narrow component is 
not significantly detected in the GSPC data, since it reduces 
X 2 by only ~ 1 for 5 additional parameters. The upper limits 
on the derived line equivalent widths are typically less than 
~ 30 eV, easily compatible with the ASCA and BBXRT 
detections. Clearly there should also be a reflection contin- 
uum accompanying this line unless a substantial portion of 
it arises from an optically thin stellar wind. We model this 
by including an additional, unsmeared, neutral reprocessed 
spectrum, assuming a fixed inclination of 60°. Again, this 
does not improve the fit (xl = 535/604), but we include 
it for completeness. Details of this fit are given in Table 1, 
where the iron abundance is frozen at its best fit value of 
1.8 x solar so that the error bars can be calculated for each 
spectrum separately. 

Several points are immediately apparent from this se- 
ries of models. Firstly, the line strength is perfectly con- 
sistent with solar, or even supersolar abundances. It is not 
anomalously weak compared to the reflected continuum as 
proposed by Ross et al (1996). The reflector is significantly 
ionized, but not to the strength that would be required for 
Auger ionization to operate (£ ~ 300 — 1000 for a power 
law of T = 1.7 and temperature of 10 6 K). The derived 
ionization state assuming an inclination of 30° is typically 
that of Fe VIII-XII, incompatible with the Auger range of 
Fe XVII - XXIII. These ionization states are obtained de- 
spite allowing for the effects of relativistic smearing. There is 
no conspiracy, whereby the high Auger ionization is masked 
by Doppler smearing and gravitational redshifts (Ross et al 
1996). Instead the line and reflection continuum are per- 
fectly consistent with each other, and with being produced 
in a disk which is weakly ionized. 
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Table 1. EXOSAT GSPC, GINGA and ASCA data from Cyg X-l. Error bars are Ax 2 = 2.7, Galactic column is fixed at 6 X 10 21 cm" 2 , 
inclination at 30° and iron abundance at 1.8x solar 



dataset 


pl r 


PL Norm 


^/ 27r rel 




Rln(Rg) 


^V^non — rcl 




GSPC 08 
GSPC 09 
GSPC 10 
GSPC 13 
GSPC 14 


i 7c; +0.02 
,o -0.01 

1.56 ±0.02 
1.58 ± 0.02 

, 7 c+0.03 
x ' ,o -0.02 

1.71 ±0.02 


1.60 
1.24 
1.16 
2.35 
1.75 


21+0- 04 
19+0- 04 

u -0.07 
H+0- 04 

0.24 ±0.05 
0.22 ±0.04 


17 +40 
1 '-16 

24^00 
30+« 

qi+60 
01 -29 

46+20 


1 '-8 

13+ 9 

1000-952 

6.0+ 6 


Q+0.07 
0+0.10 
Q+0.12 

0+o.n 

Q+0.03 


89.2/118 

92.9/98 

69.3/104 

174.8/176 

108.7/109 


GINGA 91-1 
GINGA 91-2 


1.60_ 007 
1 61+ 003 


1.61 
1.85 


0.23+0?? 
0.29«;? 7 






0+0. la 


10.2/21 
9.8/21 


GIS 1 
GIS 3 
GIS 4 
SIS 5 
GIS 6 
GIS 7 
GIS 8 
GIS 9 
GIS 10 


1.87 ±0.02 
1 84+ - 02 
1.71 ±0.02 
1.68 ± 0.03 
1.68 ±0.02 

I 44+0-03 
^ -0.02 

1.66 ±0.02 
1.59 ±0.02 
1 60+0- 03 

X ' DU -0.01 


1.79 
1.80 
1.24 
1.22 
1.89 
0.93 
1.44 
1.43 
1.62 


0.05+°g? 
06 +0 ' 03 

u - uo -0.02 

0.1 ±0.03 
09+0- 78 

u ' uy -0.04 

02+0- 09 
u - uz -o.oi 

05+0- 08 
U - UO -0.02 

06+ 06 
09+0- 09 
06+0- 11 


2300+_T 6 ° ° 
3 ± 2.4 x 10 4 

50+J°° 

3+. 32 2 x 10 4 

6 - 70 -o 9 (?i x 103 
8-4l 42 9 x 10 3 

-, 1 r'+1600 
li0 -90 
+40 

11n +2000 


16+1 
58±™ 
90i 4 20 

01 +210 
oi -26 
14+14 

q 7 +40 
°' -21 

64±J°° 


q+0.03 

0.04 ± 0.04 

0+0.05 
o 05+0-08 

n iq+0.06 

u ' lo -0.05 
q+0.03 

0.05 ± 0.05 

1O+0.05 

0.08 ±0.06 


107.8/79 
92.0/79 
142.8/79 
51.8/70 
84.4/79 
54.3/79 
108.2/79 
68.1/79 
95.8/79 



a Parameter restricted to the range 0-0.1 since it is not well constrained 



3.2 GINGA data 

The GINGA data give consistent results, though at much 
lower significance. There were 3 epochs in which Cyg X- 
1 was observed, the first of which (August 1987) was mis- 
pointed and the remaining two (in May 1990 and June 1991) 
being in an unusual mode designed to increase the spectral 
range out to high energies. In this mode the spectral resolu- 
tion is reduced by about a factor 2 since the same number 
of channels cover twice the bandpass. The resolution is then 
more like 30 per cent at iron, rather than 17 per cent in the 
usual mode. This is a crucial difference, since the broadest 
features expected from the inner disk are of order 30 per 
cent if the disk extends down to the last stable orbit, so 
these data are unlikely to be able to constrain the amount 
of relativistic smearing. Ironically, the larger bandpass is not 
actually of use either, as data beyond 30 keV are limited by 
uncertainties in background subtraction. 

We use datasets 1 and 2 (the unabsorbed spectra) 
from the June 1991 observations shown by Gierlihski et al 
(1997a). These are fit in the 4-30 keV range with the same 
model as above. The poor resolution of the data means that 
the relativistic smearing is not significantly detected, so the 
relativistic and non-relativistic reflected spectra are degen- 
erate. From the EXOSAT results and from theoretical ex- 
pectations of the amount of reflection from the companion 
star (Basko 1978) we limit the solid angle f2/27r n0 n-rci < 0.1. 
With this constraint we obtain the results shown in Table 1. 
Again the data have a line which is easily consistent with the 
amount of reflection seen. The observed ionization parame- 
ter of the reflecting material is generally low (though poorly 
constrained) even when relativistic distortions are included. 



3.3 ASCA data 

The ASCA SIS and GIS data used here are those published 
by E96, except that we neglect the October 1993 SIS spectra 
(spectrum 2 in E96) where there may be a residual problem 
with the resolution (E96). Again we use only data above 4 
keV so as to avoid the complexities of the soft X-ray excess, 
and so that the fits can be compared with the results of E96. 
These data clearly show a narrow line component as well as 
a complex edge structure which has been modelled by an 
ionized reflection (E96). However, their reflected continuum 
model did not include the self consistent iron line emission 
nor relativistic smearing. 

The SIS spectrum has the best energy resolution 
(November 1993, denoted spectrum 5 in E96) so we first ex- 
amine the line profile in these data. We fit this with the self- 
consistent iron line and continuum reflection code, together 
with a power law continuum and galactic Nh — 6 x 10 21 , 
firstly assuming an inclination of 30°. Without relativistic 
effects the fit is much worse than that shown by E96, with 
xt = 98/71, showing that the narrow line cannot be the 
corresponding fluorescence from the material giving rise to 
the observed reflection continuum, irrespective of the iron 
abundance (best fit value of unity) and ionization (£ ~ 5). 
Including relativistic effects gives a much better fit, with 
xt = 56/70, comparable with that of E96. Including a sepa- 
rate narrow line gives an even better fit, with line energy and 
equivalent width of 6.46 keV and 10 eV {xt = 48/68). How- 
ever, physically we expect a self-consistent reflected contin- 
uum to accompany this line. Replacing the line with the 
total reprocessed spectrum from neutral material where rel- 
ativistic effects are not important gives an equally good fit 
with xt = 51/69. This fit (shown in Figure 2) is very dif- 
ferent to the one derived from the EXOSAT data. The ion- 
ization is so high (5~4x 10 4 ) that the 'edge' feature seen 
in the data is actually modelled by the drop following the 
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Figure 2. The best fit to the ASCA SIS 5 spectrum assuming 
an inclination of 30°, with iron abundance free to vary between 
1 - 2x solar. The relativistic reflector has £ = 3. 0+| 7 2 X 10 4 , 
R ln = 23^ 4 J and Q/2tt = 0.09^'^, while the neutral unsmeared 
reflector has Q/2tt = 0.05^'^, for an illuminating power law 
with a = 0.68 ± 0.03 (xl = 51.8/70). 

relativistically smeared (Ri n ~ 22i? g ) ionized line, and the 
real edge is shifted outside of the observed energy range. 

The GIS spectra give similar results. We fit these using 
the most recent GIS response matrices (as of March 1995, 
gis2v4-0.rmf) as opposed to the earlier response used by 
E96, and results are detailed in Table 1. Again, the derived 
ionization paramter is generaly very high - only the GIS 4 
and GIS 9 spectra give results that are anything like the low 
ionization EXOSAT and GINGA solutions. 

3.4 Inclination 

This general mismatch between the ionization of the repro- 
cessor derived from ASCA and EXOSAT/ GINGA data is 
suggestive of a systematic problem in the spectral model 
fitting, most likely due to the difference in bandpass. In EX- 
OSAT and GINGA the high energy continuum shape of the 
reprocessed spectrum helps constrain its ionization, whereas 
in ASCA only the iron features can be used. The detailed 
shape of the line and edge is a strong function of inclination 
as well as of ionization. At high inclinations Doppler shifts 
prevail over gravitational and transverse redshift, shifting 
the line and edge to higher energies as well as giving sub- 
stantial broadening. This is to zeroth order the same effect as 
ionization. However, the detailed shape of a low inclination, 
ionized reflection spectrum is rather different to a less ion- 
ized reflection spectrum at higher inclination, so these two 
parameters can be disentangled given good enough data. 

Figure 3a shows how x 2 varies for simultaneous spectral 
fitting of the five EXOSAT datasets as the disk inclination 
is changed from 30 — 66° . Clearly there is a significant pref- 
erence in the data for an inclination higher than 45° (cosine 
smaller than 0.7). The 3<r (A^ 2 = 9) limit to the inclination 
from this modelling is i > 40°. This is rather higher than 
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Figure 3. The change in x 2 (goodness of fit parameter) as a 
function of inclination of the reflecting material. Panel (a) shows 
results from a simultaneous fit of all the EXOSAT GSPC spec- 
tra, with the solid line showing fits to the data from 4-20 keV 
while the dotted line uses the full 2-20 keV bandpass including a 
(diskblackbody) model for the soft excess. The latter have 72 sub- 
tracted from them in order to appear on the same scale. Panel (b) 
shows the ASCA GIS 3 (solid line) and ASCA SIS 5 (dotted line) 
4-10 keV spectral results, where the latter have been boosted by 
20 in order to appear on the same scale. Clearly the data which 
have significant bandpass beyond 9 keV favour inclinations higher 
than 30° . 



the most probable inclination inferred from optical studies 
of 28° - 38° (Gies & Bolton 1986), but within the limit set 
by the most recent work of < 55° (Sowers et al 1998). At 
the best fit GSPC inclination of 53° the iron abundance is 
2±°- b 4 x solar (xl = 522/604). We check that this is not an 
artifact of ignoring the soft X-ray excess by including data 
down to 2 keV and approximating it with a diskblackbody 
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spectrum (XSPEC model diskbb). The dotted line in figure 
3a shows that the data still significantly favour the higher 
inclination solutions. 

Despite their higher spectral resolution, the ASCA SIS 
are completely unable to constrain the inclination (Figure 
3b, dotted line). Figure 4a shows this degeneracy of SIS 
solutions in more detail. At low inclinations the data select 
uniquely the high ionization solution (shown in Figure 2), 
while at high inclinations an equally good solution can be 
found at both high and low ionization. Figures 5a and b 
show these high inclination solutions - the low ionization 
one has the 'pseudo-edge' in the data matched by the real 
edge while at high ionization it is matched by the drop from 
the broadened ionized line. Figures 5a and b also show the 
extension of these models out to 20 keV. It is obvious that 
these solutions can be distinguished with higher energy data. 

The GIS data have just enough extra bandpass at high 
energies (above 9 keV) to break this degeneracy. The dashed 
line in figure 4b shows the change in \ with ionization for 
cos i — 0.86 and 0.6 for the GIS spectrum 3 of E96. Not only 
is the ionization uniquely determined for a given inclination, 
the data also are able to distinguish between the high ion- 
ization/low inclination and low ionization/high inclination 
solutions, clearly preferring the latter. The correspondence 
with the GSPC data is clear, so we fix the inclination of 
the relativistic reflector at 53° (cosine of 0.6) and re-fit all 
the data, using the best fit GSPC iron abundance of twice 
solar. These results are shown in Table 2. For the EXOSAT 
GSPC (and GINGA) data the fits are only subtly different - 
the derived ionization state is somewhat lower while the in- 
ner radius is somewhat larger, both because of the stronger 
doppler effects in the higher inclination models. However, 
for the ASCA data the high inclination solutions are very 
different from those in Table 1. Firstly, they are statistically 
better - the total \ 2 from all the GIS data for the high incli- 
nation solution is 698.5/632 as opposed to 753.4/632 for the 
low inclination one. Secondly, the derived ionization drops 
dramatically, and now corresponds to that seen in the EX- 
OSAT GSPC data. All the data now give similar derived 
parameters, irrespective of bandpass. Only the SIS 5 and 
GIS 6 spectra allow a high ionization solution. 

Unlike the variable gain EXOSAT GSPC data, the GIS 
data are all taken in the same mode, so their spectral residu- 
als can be co-added in order to illustrate the features which 
make one model better than another. The data are fit si- 
multaneously, with only the reflection parameters of iron 
abundance (constrained to be between 1-2 x solar) and in- 
clination fixed between the datasets. The top panel in Figure 
6 shows the results from fitting a reflector which is ionized 
but not relativistically smeared for an inclination of 30°. 
There are clear broad residuals present around the iron line 
energy. A much better fit is obtained if the 30° ionized re- 
flector is allowed to be relativistically smeared (and a non- 
smeared, unionized reflector is also added), as shown in the 
middle panel. However, the bluest features of the residuals 
are not able to be adequately fit in these models, despite the 
rather high derived ionizations (see table 1), due to the pre- 
dominant redshift at 30°. Higher inclination solutions give 
a stronger doppler shift to the line, so matching the data 
rather better (bottom panel). 

For the high inclination solutions, all the data show 
that the covering fraction of the relativistic reflector is al- 
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Figure 4. The change in the goodness of fit parameter \ 2 as a 
function of ionization state of the reflector for inclinations of 30° 
and 53° . Panel (a) shows the results for the SIS 5 spectrum. For 
an inclination of 30° the data uniquely select the high ionization 
solution while at 53° there are equally good solutions at both low 
and high ionization. All these three solutions are of comparable 
statistical likelihood. Panel (b) shows the contrasting situation 
derived from the GIS 3 spectrum. Here the data uniquely select 
the high inclination, low ionization solution. 



ways much less than unity, the derived inner disk radius is 
generally significantly larger than 6i? g , and its ionization 
state is rather low. The solid angle subtended by the narrow 
reprocessor is always consistent with being the £l/2iv < 0.1 
expected from the companion star. This should vary with or- 
bital phase (Basko 1978), although the data here give error 
bars on the derived parameters which are too large to con- 
strain this. Orbital variability in the narrow reflected com- 
ponent would provide an explanation for the observed 4% 
peak to peak quasi-sinusoidal variability seen in the high en- 
ergy BATSE (20-100 keV) lightcurve of Cyg X-l (Paciesas 
et al 1997). The reflected component should be small (basi- 



© 0000 RAS, MNRAS 000, 000-000 



Relativistic Distortions in the X-ray Spectrum of Cyg X-l 7 




(a) 



> 

34 



S3" 
(X o 



^ 2 



5 10 
channel energy (keV) 



20 




(b) 



34 ^ 
m O 



5? 
Ph o 



^2 5 10 20 

channel energy (keV) 

Figure 5. Two equally acceptable decompositions of the SIS 
5 spectrum of Cyg X-l, assuming an inclination of 53° with 
iron abundance free to vary. Panel (a) shows the low ioniza- 
tion reflected continuum, with a relativistic disk of ionization 
£ = 0.04lf, 5 04 , R in = 14 ± 4 and Q/2tv = 0.45±°;*° {xl = 
49.7/69), for an illuminating power law with a = 0.731^ and 
neutral unsmeared reflection Q/2w = 0.05~^q'q®. The iron abun- 
dance is completely unconstrained between 0.5 — 3x solar. Panel 
(b) shows the high ionization solution, where the relativistic disk 

nth 



+4.0 



has £ = 3+2 6 X 10 4 , R in = 55Zlo and Q/2w = 0.18 
Q/2-7T = 0.07 ± 0.06 of unsmeared neutral reflection and a power 
law with a = 0.70 ± 0.04 {xl = 50.0/69). 



cally zero) at phase when the companion star is in front of 
the black hole, and have its maximum of £l/2n ~ 0.1 viewed 
at an inclination of 90 — i at phase 0.5, giving a 5% peak-to- 
peak variability in the 20-100 keV flux. We do not require 
an additional narrow reflected component from the outer 
disk, contrary to the conclusions of E96. Their fits to these 
data gave some spectra which required a much larger narrow 
component than could easily be explained by the compan- 
ion star alone (e.g. GIS data from November 1994, spectrum 
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Figure 6. Co— added \ 2 residuals from fitting the 8 GIS spectra. 
The top panel shows results from fitting non— smeared reflection 
to the data. There are clearly broad residuals left around the iron 
line energy. The mid panel shows that these are much reduced 
by allowing the reflector to be a relativistic disk (including non- 
smcarcd, unionized reflection from the companion star). However, 
for the assumed inclination of 30°, the highest energy residuals 
are not matched by the models. The lower panel shows that these 
can be fit at higher inclinations, where doppler blueshifts can 
become stronger than the prevailing rcdshifts. The feature at 4.7 
keV in all the plots is a residual from the changing response of 
the GIS detectors at the Xe L i 



10), even allowing for twice solar iron abundances. The rea- 
son for the difference is that the fit is somewhat sensitive 
to changes in the model, and we note that our description 
gives better \ 2 (by up to 20) than the fits of E96. 



4 IONIZATION AS A FUNCTION OF RADIUS 

The data clearly prefer a solution in which the covering frac- 
tion of the reflector is extremely low, typically only 20-30 
per cent of what would have been expected from a disk which 
covered half the sky as seen from the X-ray source, and the 
derived inner radius is generally inconsistent with the last 
stable orbit. Both these results can easily be explained in 
a scenario where the optically thick disk terminates before 
the last stable orbit, with the inner radii instead being filled 
with X-ray hot optically thin(ish) plasma. However, a series 
of papers (e.g. Ross et al 1996) have explored the possibility 
that high ionization of the reflector supresses the derived 
amount of reflection. We have shown above that this does 
not work for Cyg X-l when the reflector is modelled by 
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Table 2. EXOSAT GSPC, GINGA and ASCA data from Cyg X-l. Error bars are Ax 2 = 2.7, Galactic column is fixed at 6 X 10 21 cm" 2 , 
inclination is 53° and the iron abundance is 2x solar 



dataset 


pl r 


PL Norm 








^/ 2n non — rcl 




GSPC 08 
GSPC 09 
GSPC 10 
GSPC 13 
GSPC 14 


1.77 ±0.02 

1 cro+0.02 
, co+0.02 

i-os_o.cn 

1 yo + 0.02 

1 74+0.01 
' -0.02 


1.65 
1.28 
1.16 
2.46 
1.85 


n 34+0.05 

26+ - 09 
U - ZD -0.12 
13 +0.05 

34+ - 12 
0.39 ±0.06 


2+ 17 
Q+0.3 

21+ 70 
0+3 
Q+0.4 


28t 22 

1000-970 
1411 


Q+0.14 

05+ 011 
0+0-13 

09+ 10 
Q+0.05 


86.1/118 
96.3/98 
69.7/104 
172.6/176 
98.0/109 


GINGA 91-1 
GINGA 91-2 


1.63 ± 0.03 
1.63 ± 0.03 


1.60 
1.92 


no + 0.15 

44+0-10 
u -^-0.21 


4lf 

o+45 
a -3 


4oir 

30lS 


0.1-0.10 

0+0. la 


12.1/21 
8.8/21 


GIS 1 
GIS 3 
GIS 4 
SIS 5 
GIS 6 
GIS 7 
GIS 8 
GIS 9 
GIS 10 


1.91 ± 0.02 

1.91 ± 0.02 

1.74 ±0.02 

1.74 ±0.04 

1 69+ 03 
x - oy -0.02 

1.50 ±0.03 
1 66+ 02 
1.58 ±0.02 
1.61 ± 0.02 


1.89 
2.08 
1.28 
1.46 
1.94 
1.01 
1.45 
1.42 
1.63 


0.36+" 2 

— 0.13 
^o+O.OO 
u - oa -0 .17 

0.31 ±0.08 
54+ - 17 

o.io±°:SS 

0.481HI 

09+ - 05 

04+ - 13 

04+ - 09 
u - u ^-0.02 


2.8+f- 

— 2.7 

02+ 17 
0+ 1 

0+20 

n n +le5 
JU -30 

0+io 

50 +2000 

0+oo 
700 +1700 


20 -t 

9 9 +6 
a -1.9 

6.0+ 2 
13+* 

i4±i 

60±i™ 

iool 220 


°- 01+ ooi 

aoslgiS 

0.09 ±0.03 
0.05 ±0.03 
0.12 ±0.03 

Q+0.02 

05+ 005 
U - UO -0.04 

12+ - 04 
u ' lz -0.12 

09+ 04 
u - uy -0.05 


93.5/79 
69.8/79 
134.4/79 
49.9/70 
81.4/79 
45.7/79 
108.4/79 
71.7/79 
93.6/79 



a Parameter restricted to the range 0-0.1 since it is not well constrained 



Table 3. GSPC 08, SIS 5 and GIS 3 spectra fit to the radial ionization model with £(r) oc r~ 1,5 . Error bars are A\ 2 = 2.7, Galactic 
column is fixed at 6 X 10 21 cm -2 . Inclination is fixed at 30° and 53° for each spectrum, while the iron abundance is allowed to vary between 
1-2 X solar (so there is an extra degree of freedom compared to table 2). 



Spectrum PL T A Fc Cos(i) U/2-K iei £(-Rm) Rin ^/27r non _ rcl xl 



GSPC 08 


i 77 +0.01 
x - ' 1 -0.03 


2.0 


0.60 


34+ - 04 
U " ;i4 -0.14 


n + 50 
Z -2 


27 +i7 


i+o.i2 

u ' ul -0.01 


86.2/117 


GSPC 08 


1.77 ±0.04 


1.2 


0.86 


0.221H 4 


fiC ; + 160 


201 7 3 


0+0.09 


88.5/117 


GIS 3 


I ni+0.09 
1 - a -0.02 


2 


0.60 


o.53t-2:« 


0-Ollaoi 


9-5t 3 2.5 


09+ 07 
u - ua -0.04 


69.7/78 


GIS 3 


1 85+ 001 
1 ' oo -0.02 


2 


0.86 


0.06 ±0.01 


3-0+1.1 x 104 


o 4 +27 
° 4 -12 


0+0.06 


89.7/78 


SIS 5 


i 7C :+0.13 
1 ''°-0.04 


2 


0.6 


5fi+ - 59 

U.OD_ 12 


1+ 30 


13+6 


06+ 08 
u ' UD -0.03 


49.4/69 


SIS 5 


1 69+ 05 

-0.03 


2 


0.6 




l.40l 2 ;3 7 x 106 


17 34 
1 '-10 


0-Htaoe 


50.7/69 


SIS 5 


1 67+ 05 
J -°'-0.03 


2 


0.86 


1 1 +°° 


1.8lf 5 x 10 5 


14+16 


0-05lS° 04 


52.0/69 



material with constant ionization parameter. However, per- 
haps there is a radial range of ionization states. If there is 
an extremely ionized inner disk then it produces no spectral 
features, so neither its covering fraction nor its relativistic 
smearing can be seen. 

We test this idea by dividing the disk into 10 radial 
zones, again with the illumination oc r -3 , and the ionization 
f(r) oc r -1,5 as appropriate for a gas pressure dominated 
disk. Due to computational time we only demonstrate this 
model on a single spectrum from each instrument: EXOSAT 
GSPC 08, ASCA GIS 3 and ASCA SIS 5. The iron abun- 
dance is a free parameter (within the range 1 — 2x solar), 
and we include the possibility of a neutral reflected compo- 
nent from the star/outer disk. The fits are repeated for fixed 
inclination of cosi = 0.86 and 0.6. 

Details are given in Table 3, and are not significantly 
different from those derived from the single ionization model. 
The inferred covering fraction for the disk is still low, and the 
disk inner radius is inconsistent with the last stable orbit. 
We explicitally searched for a high ionization solution for the 
cos i — 0.6 fits, but these always gave a significantly worse 
X 2 except for the SIS 5 data. 



The reason for the failure of the multi-zone ionization 
model is that the mean ionization state observed is rather 
low. Between this and the extreme ionization states required 
to make the reflector invisible are the intermediate/high ion- 
ization states. These give a strong and distinctive spectral 
signature, especially those where iron has only one or two 
electrons left (H and He-like, FeXXVI and XXV) which pro- 
duce intense line emission at 7.0 and 6.7 keV, and (more im- 
portantly) a deep edge at 9.2 and 8.6 keV (Ross & Fabian 
1993). Even allowing for strong distortion by relativistic 
smearing (Matt et al., 1993a; 1996; Ross et al., 1996) these 
features are not seen in the data. It is not then possible to 
go smoothly from the observed rather low ionization states 
to the extreme ionizations required to hide reflection. Only 
if there is an abrupt ionization transition can this model be 
sustained. 

Again, we note that good data at energies > 9 keV are 
required in order to constrain the amount of highly ionized 
reflecting material. Radial ionization models with high inner 
disk ionization can be fit to the SIS data for any inclination, 
but are again ruled out by the GIS 3 spectrum (see table 3). 
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5 DISCUSSION 

From the above fits it is clear that the relativistic smear- 
ing of the reflected spectrum is significantly detected in Cyg 
X-l, confirming the presence of an optically thick material 
at small radii - the putative accretion disk. The covering 
fraction of the disk is substantially less than unity (see Ta- 
bles 1-3), irrespective of inclination, ionization and radial 
gradients in ionization. 

The ionization of the reflector and its inclination to the 
line of sight can be constrained with moderate spectral reso- 
lution data and bandpass extending beyond 9 keV. The data 
show strong preference for a low ionization reflector inclined 
at rather more than 30° to the line of sight: the ASCA GIS 
3 spectrum gives a range of 41° — 61°, while the EXOSAT 
GSPC data imply angles greater than ~ 47° (both at the 
90% confidence limit). These inclinations are determined by 
the detailed shape of the iron line in the data: for low in- 
clinations the theoretical line profile is strongly affected by 
gravitational redshift so that higher ionization solutions are 
needed to fit the data. But in higher ionization solutions 
there is a large energy separation between the line and edge 
(however smeared they become) which is not seen in the 
data. The spectral features observed require that the line 
and edge blend together, i.e. that the ionization is moder- 
ate to low. The reflector then is required to be viewed at a 
moderate inclination angle so that the doppler blueshift can 
counteract the gravitational redshift. However, these inclina- 
tion are rather larger than the often quoted probable range 
of 28° — 38° (Gies & Bolton 1986) derived from optical data, 
although within the firm upper limit to the inclination from 
these studies of i < 55° (Sowers et al 1998). The lack of a well 
determined optical orbital solution is due to the complexity 
of the accretion flow from the supergiant companion in this 
system. The companion probably does not completely fill 
its Roche lobe, but instead has a strong stellar wind which 
is gravitationally focussed by the Roche potential (see e.g. 
numerical simulations by Blondin, Stevens & Kallman 1991) 

While the iron line is strongly distorted by the relativis- 
tic effects, the extremely broad components from material 
at 6i?g are not present at the level expected from a disk 
extending down to the last stable orbit with line emissivity 
oc r~ 3 . Such a model is only possible if the inner disk has 
an abrupt increase in ionization state to the level at which it 
becomes completely ionized so produces no atomic spectral 
features. We stress that a continuous (power law) ionization 
gradient as expected from the intense X-ray illumination is 
insufficient to do this, since between the observed rather low 
ionization state of the reflector and the extreme ionizations 
required to make the disk completely reflective there are the 
high ionization states, where iron has only 1 or 2 electrons 
left. These ions produce copious Ka line emission at 6.7 and 
7.0 keV, and strong edges at 8.8 and 9.2 keV, which are not 
seen in the data. Similarly, these high ionization states are 
precisely the ones expected from collisional ionization from 
the accretion disk temperature. The observed soft excess is 
dominated by a component at 0.1 — 0.2 keV which cannot 
completely ionize iron. Even the secondary soft emission at 
kT ~ 1 keV is far removed from the inner K shell energy 
of 8.8 and 9.2 keV for He and H like iron. This then rules 
out the idea of a simple correspondence between AGN and 
GBHC where both have disks which extend down to 6R R 



but in the GBHC the disk is ionized simply because of the 
higher accretion disk temperature. 

The optically thick disk must do something at the de- 
rived 7?i n > 67? g (assuming that the line emissivity is oc r~ 3 ; 
Zycki ct al 1998b). Either the inner disk disappears com- 
pletely or it jumps dis continuously to a completely ion- 
ized state at this radius. Similar results are seen for the 
low/hard state GBHC transient sources (Zycki et al 1997; 
1998ab), showing that this is a robust feature of low/hard 
state GBHC spectra. 

Figure 7 shows the plot of spectral index versus the de- 
rived covering fraction, ionization and inner radius of the 
reflector for the reflector models inclined at 53° (not includ- 
ing the poorly constrained GINGA data). The covering frac- 
tion is not consistent with a constant (\l = 60/14), showing 
that the geometry must be changing. Apart from the hardest 
spectral index data (GIS 7, which incidentally is the spec- 
trum taken furthest off axis, and so has a large vignetting 
correction, E96) the data are consistent with a trend in 
which a hard spectral index is associated with smaller re- 
flected fraction, similar to results from some of the X-ray 
transients (Zycki et al 1997; 1998ab). This can be qualita- 
tively explained by a geometry in which there is a variable 
radius inner hole in the cool disk which is filled with X-ray 
hot plasma, although the error bars on the radius derived 
from relativistic smearing are too large to be able to con- 
strain these ideas. A disk which extends further down into 
the gravitational potential could subtend a larger solid an- 
gle to the source. The soft photon flux from the disk would 
also increase, both because the disk luminosity is dependent 
on its inner radius, and also from the larger fraction of hard 
X-rays which can be reprocesssed. This would give stronger 
Compton cooling of the hot electrons and so a steeper X-ray 
spectrum. 

The typical derived inner radius of the disk is Ri n — 
20i?g. The luminosity released by viscous dissipation from 
20i?g to 6i?g is rather less than 0.75 that released from oo to 
20i?g (see e.g. equation 5.19 in Frank, King & Raine 1992). 
Simple energetics then give the expectation that the X-ray 
hot plasma should carry ~ 0.75 x the luminosity of the soft 
photons from the disk. However, the opposite is observed: 
the hard component (the X-ray hot plasma) is over twice as 
luminous as the soft component from the disk e.g. di Salvo 
et al., (1998). 

However, this geometrical interpretation is clearly not 
unique, and many others have been proposed (see e.g. 
the review by Poutanen 1998). If the hot X-ray emitting 
plasma forms a (plane parallel) corona above the disk (e.g. 
Haardt & Maraschi 1993) then the reflected spectrum has to 
cross the corona before escaping to the observer. A fraction 
cxp(— t / cos i) (where r is the optical depth and i the inclina- 
tion of the disk-corona to the observer) will interact again 
with the hot electrons in the plasma. The reflected spec- 
trum then forms another source of seed photons for Comp- 
ton scattering, and the scattered fraction joins the Compton 
continuum and so is unrecognisable as reflected flux. Such 
models were shown to work for Cyg X-l data up to 20 keV 
(Haardt et al 1994), but they predict that there should be a 
strong break in the spectrum beyond this energy. Seed pho- 
tons from the disk are not isotropic, so the first few Compton 
scattering retain the imprint of this anisotropy and give a 
somewhat flatter spectrum which then joins onto the more 
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Figure 7. The relativistic reflection spectrum parameters as a 
function of the intrinsic spectral index for the model with reflector 
inclination of 53° (table 2). The top panel shows the amount of 
reflection, while the second panel shows its ionization state and 
the third shows its inner radius. 

isotropic higher order scatterings (Haardt & Maraschi 1993). 
This anisotropy break is not present in broad band spectra 
from Cyg X— 1, where the data extend past 20 keV, which 
rules out such models (Gierlinski et al 1997a). 

The lack of an anisotropy break seems to be an in- 
surmountable problem for any model which has the X-ray 
emission region above a disk, including the 'active regions' 
models which may be more appropriate if the X-ray regions 
are energised by localised magnetic reconnection events (e.g. 
Haardt, Maraschi & Ghisellini 1994). However, these models 
are very attractive given that accretion disk viscosity is now 
thought to be driven by magnetic reconnection. But even if 
the anisotropy break can be masked (perhaps by the elec- 
trons also being anisotropic or by a distribution of electron 
temperatures) there are further problems for models where 
the X-ray region is above the disk. If edge effects are not im- 
portant (approximately plane parallel geometry) then about 
half of the hard X-rays are intercepted by the disk. The ma- 
jority (whatever is not reflected) of this flux is thermalised in 
the disk and re-emerges as soft photon seeds for the Comp- 
ton scattering. These soft photons all have to pass through 
the corona in a plane parallel geometry, resulting in strong 
cooling. Without a strong anisotropy break then the pre- 
dicted spectra are too soft to match the observed low/hard 
spectra from Cyg X-l (Dove et al 1997, Poutanen et al., 
1997). This forces us to the 'active regions' models, where 
not all flux of soft photons from the disk will pass through 
the corona as edge effects are important. But in these models 
there is then no overlying corona that can reduce the amount 
of reflection seen. If edge effects are important, then much 
of the reflection is produced from regions of the disk that 
are not underneath the corona. Reflection is not then su- 



prcsssed by Compton scattering in the corona, and there is 
the strong prediction that we should see 57/27r ro i = 1, which 
is plainly not the case. 

These are serious problems to be faced by all models 
which have the X-ray region above the disk. The continu- 
ous disk — corona models can reduce the amount of reflection 
seen, but predict a large anisotropy break which is strongly 
ruled out by the data. Even if there is some way to mask the 
anisotropy break then these models produce spectra which 
are too soft to match what is seen. The 'active regions' mod- 
els can produce hard spectra, but cannot supress reflection 
and again predict an anisotropy break which is not seen. 
Compared to these drawbacks, the central source geometry 
seems an attractive option. 



6 EQUILIBRIUM ACCRETION FLOW 
SOLUTIONS 

There are several solutions of the equations for accretion of 
material onto a compact object. The most well known of 
these is that of Shakura & Sunyaev (1973; hereafter SS), 
which is the geometrically thin, optically thick accretion 
disk. In this solution, the gravitational energy released is 
radiated away locally as (quasi)blackbody radiation. How- 
ever, other solutions are also possible: the energy can instead 
be radiated by Compton scattering in an optically thin(ish), 
geometrically thick(ish) flow (Shapiro, Lightman & Eardley 
1976). Alternatively, the energy released need not be effi- 
ciently radiated. It can be carried along with the flow (ad- 
vected) and swept into the black hole. These inefficient, ad- 
vectively cooled flows can either be optically thin (ADAFs: 
Narayan &: Yi 1995; Abramowicz et al 1995), or optically 
thick (Abramowicz et al 1988), depending on the mass ac- 
cretion rate. In general the solutions coexist: at any radius 
there are several different equilibrium flow configurations for 
a given mass accretion rate and viscosity (see e.g. Chen et 
al., 1995). 

6.1 Shakura— Sunyaev Accretion 

Why should the inner disk radius change, and not extend 
down to 6-Rg ? The 'standard' accretion disk model devel- 
oped by SS gives the flow structure in the limit when the 
gravitational energy is released in optically thick material. 
The cooling is very efficient (blackbody radiation), so all 
the energy released locally can be radiated locally, giving a 
cool, geometrically thin disk structure. This becomes unsta- 
ble when the total pressure is dominated by radiation pres- 
sure (Shakura & Sunyaev 1976), giving a natural inner disk 
cutoff. For steady state models these predict that the inner 
radius should move outwards as the luminosity increases (as 
the radiation pressure dominated region increases: SS), in 
conflict with the observed trend in the transient GBHC to 
show higher reflected fractions and smaller inner disk radii 
at higher luminosities (Zycki et al 1997; 1998a,b). 

Another, more general problem with the SS accretion 
disk structures is that such models give temperatures of or- 
der 1 keV for GBHC at high accretion rates but are quite 
unable to explain the observed (hard or soft) power law 
tail out to energies > 100 keV. Either there are parts of 
the disk flow in a different configuration to that of SS, or 
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there is some non-disk structure such as a corona powered 
by magnetic reconnection (e.g. Haardt, Maraschi & Ghis- 
ellini 1994). However, this latter possibility (a cool disk un- 
derneath a plane parallel or hemispherical X-ray emitting 
corona) seems to be ruled out by detailed spectral fitting 
(see Section 6). Instead it seems that the spectrum implies 
a geometrical constraint that the X-ray emitting plasma and 
disk occupy more or less separate regions, such as an outer 
disk and inner (quasi-spherical ?) X-ray source (Gierlihski 
et al 1997a, Dove et al 1997, Poutanen et al., 1997). 

6.2 Advective Accretion 

Recently there has been much excitement about the possi- 
bility that another stable solution of the accretion flow may 
explain the hard X-ray data. Below a critical accretion rate, 
rh < rh cr it, a stable, hot, optically thin, geometrically thick 
solution can be found if radial energy transport (advection) 
is included (see e.g. the review by Narayan 1997). The key 
assumption of these models is that the protons gain the en- 
ergy from viscous processes, while the electrons only acquire 
energy through interactions (electron-ion coupling) with the 
protons. This coupling timescale can be rather slow com- 
pared to the accretion timescale, so protons can be accreted 
into the black hole, taking some of the energy with them. 
The energy that the electrons do manage to obtain is radi- 
ated away via cyclotron/synchrotron emission (on an inter- 
nally generated magnetic field), bremsstrahlung, or Comp- 
ton scattering of the resultant spectra of these two processes. 
In contrast to the SS accretion flow models, there is no cold 
disk in the inner regions, so no strong source of soft seed 
photons for Compton scattering, hence the resulting X-ray 
spectra are typically hard. 

Such flows were proposed to explain the hard and very 
faint X-ray spectra seen from the transient X-ray GBHC in 
quiescence (e.g. Narayan, McClintock & Yi 1997), and then 
extended by Esin, McClintock & Narayan (1997) to cover 
the whole range of luminosity seen in the GBHC transients. 
As the mass accretion rate increases from rh « rh cr i t to 
rh ~ m cr it the flow density increases, so the electron-ion cou- 
pling becomes more efficient and the fraction of energy the 
electrons can drain from the protons increases. This increase 
in radiative efficiency continues to rh — m cr it, where only 
~ 35 % of the heat energy is advected. At higher rh > m cr i t 
the advective flow collapses into an SS disk. Esin et al (1997) 
identify this change from a hot advective flow to a cool SS 
disk as the origin of the hard / soft spectral transition seen 
in GBHC, and some observational support for this is indeed 
seen (Zycki et al., 1998a, Cui et al., 1998, Gierlihski et al., 
1997b). 

These models have been specifically applied to Cyg X-l 
by Esin et al (1998). Their calculated spectra approximately 
match the overall 2-500 keV GINGA-OSSE data from Cyg 
X-l for Ri n ~ 200i?g, but uncertainties in their model may 
extend this to Ri n > 40i? g . This is still larger than most of 
the constraints from relativistic smearing (table 2). Our data 
indicate that the optically thick material extends much fur- 
ther into the gravitational potential than anticipated by the 
advective flows. Similar results are found for the comparison 
of the models of Esin et al (1997) against the observed hard 
state spectra from the GBHC transient source Nova Muscae 
(Zycki et al., 1998a). 



6.3 Shapiro— Lightman— Eardley flows 

An alternative X-ray hot solution to the cool SS accre- 
tion flow was discovered by Shapiro et al., (1976, hereafter 
SLE). These solutions are actually very closely related to 
the ADAFs discussed above in that they assume that the 
protons gain most of the energy from viscous processes, and 
then heat the electrons by electron-ion coupling, resulting in 
an optically thin(ish), geometrically thick(ish) flow. This so- 
lution of the accretion flow equations merges into the ADAF 
branch at a critical mass accretion rate (see Chen et al 1995; 
Esin et al 1997; Zdziarski 1998). The key difference is that 
on the SLE branch, Compton cooling dominates over advec- 
tion, while the opposite is true on the ADAF branch. There 
is no known possible hot solution for mass accretion rates 
above the critical one at which the ADAF and SLE solutions 
merge. 

In its original form, the SLE solutions were shown to 
be viscously stable but thermally unstable (Pringle 1976). 
An increase in the heating rate is not balanced by an in- 
crease in the cooling because the electron-ion collisons be- 
come less efficient as the temperature increases. The key 
to the growth of interest in ADAFs is that they are sta- 
bilsed against this due to the dominance of advective cool- 
ing. However, the SLE is not necessarily thermally unstable. 
The stability analyses which have been done generally as- 
sumed that the viscous timescale is long compared to the 
thermal timescale, so that the density of the region can- 
not change during the perturbation. This is not true for a 
geometrically thick(ish) flow. Many stability analyses also 
neglected advective cooling, and/or external Compton cool- 
ing on seed photons which vary. We urge further study of 
the thermal stability of these flows. 

Understanding the stability (and hence the existence) of 
these flows is important since the SLE composite geometry 
is rather similar to that inferred from the data, with a hot, 
quasi-spherical plasma cooled by radiation from an external 
SS accretion disk. The advantage the SLE solution has over 
the ADAF is that it radiates efficiently, which is generally re- 
quired from the rapid observed variability of accreting black 
holes. However, the disk transition again occurs where radi- 
ation pressure dominates over gas pressure, so these models 
again apparently predict that the inner disk radius is larger 
at high mass accretion rates, opposite to the observed trend. 



7 UNSTABLE ACCRETION IN AN SS DISK 

The thermal instability above is produced wherever an in- 
crease in the mass accretion rate is not balanced by a cor- 
responding change in the cooling rate. Similarly, a disk is 
viscously unstable if a small increase in the external mass ac- 
cretion rate cannot be balanced by the viscosity to produce a 
corresponding change in the mass accretion rate through the 
disk. These two are related in that flows which are viscously 
unstable are also necessarily thermally unstable. 

The ionization of Hydrogen produces a disk transition 
which is viscously (and thermally) unstable. Below temper- 
atures of ~ 5000 K the opacity is dominated by molecular 
lines, while above this temperature there is a marked in- 
crease in opacity due to bound-free ionization of Hydrogen. 
However, for temperatures beyond ~ 10 4 K the bound-free 
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opacity drops since hydrogen is then mostly ionized. This 
gives a local maximum in opacity at ~ 6000 — 10 4 K (e.g. 
Hure et al 1994). In this range an increase in rh gives rise to 
a large increase in opacity, trapping the radiation and rais- 
ing the central (equatorial plane) temperature of the disk. 
This in turn raises the sound speed, c s , and the kinematic 
viscosity, v — ac s H, leading to an increase in the rate at 
which the disk material is accreted onto the compact object 
(see e.g. the review by Osaki 1996). 

If this instability can propagate throughout the disk 
then this leads to two distinct states. In quiescence the disk 
is cool and has low kinematic viscosity - generally so low 
that the rate at which material is added to the disk from 
the external accretion source is faster than the rate at which 
mass can be viscously transported inwards from the outer 
disk. Mass builds up at the outer edge of the disk, and the 
structure is not that of a cool steady SS disk. Eventually, the 
outer disk mass is so large that hydrogen can be ionized, and 
then the disk switches to a hot state where the kinematic 
viscosity is high and steady state can be achieved. 

This instability is thought to be the trigger for the out- 
burst in the transient black hole systems (see e.g. Lasota, 
Narayan &: Yi 1996), though it was first discovered in the 
context of the quiescent/out burst behaviour of disk accret- 
ing white dwarfs (see e.g. the review by Osaki 1996). How- 
ever, one major problem with these models is that the ac- 
cretion rate inferred in quiescence (as evinced by the X- 
ray emission) is rather larger than the maximum accretion 
rate for which the cool quiescent disk can be sustained if 
the disk extends all the way down to the compact object 
(Lasota et al., 1996). This led to the idea of a 'hole' in 
the inner disk in quiescence, perhaps formed by evaporation 
(Meyer & Meyer-Hofmeister 1994). However, another possi- 
bility which has yet to be explored in this context is that the 
inner disk in quiescence is simply optically thin. The surface 
density profile for the non-stationary cool quiescent disk is 
peaked to the outside of the disk, so it is very easy for the 
inner regions to be optically thin (e.g. Cannizzo 1998, figure 
7), especially as the absorption opacity of the cool disk is low 
so even if the photon is scattered many times, there is still 
a high probability that it can ultimately escape (effectively 
optically thin). This means that the inner disk is unlikely 
to be able to thermalise the energy released locally, and so 
may not radiate (i.e. cool) efficiently. The material would 
then heat up, and turn into something like an advective or 
externally cooled SLE flow which would appear to emanate 
from a (non stationary) optically thick disk. 

However, the main problem with such models is that 
while this may describe the situation in quiescence, it seems 
highly unlikely that it has any relevance to the high mass 
accretion rates required for Cyg X-l and the other X-ray 
bright transient GBHC. Another instability that can give 
the same effect but at high mass accretion rates is required. 
One obvious possibility is the instability which occurs when 
radiation pressure starts to dominate in the disk. The in- 
ternal pressure (and hence the sound speed and kinematic 
viscosity) then rises very rapidly (oc T 4 ), and the mass ac- 
cretion rate through the radiation pressure dominated part 
of the disk can then increase faster than the increase in 
external mass accretion which triggered the rise in temper- 
ature, so leading to depletion of the inner disk. At higher 
mass accretion rates the flow is again stabilised though the 



excess radiation being carried along with the flow and ad- 
vected into the black hole (Abramowicz et al 1988). Models 
of such flows are beginning to be developed (Szuszkiewicz & 
Miller 1998) in response to observations of rapid inner disk 
instabilities in the superluminal accreting black hole candi- 
date GRS1915+105 (Belloni et al 1997), but are currently 
in their infancy. It remains to be seen whether models of a 
non-stationary flow inside the radiation pressure dominated 
regime could reproduce what is seen in Cyg X-l and other 
low state GBHC spectra. 



8 INHOMOGENEOUS ACCRETION 

It seems highly unlikely that accretion flow can change 
smoothly between an optically thick cool solution and an 
optically thin hot solution. Some transition region where 
cool clumps are embedded in hot plasma seems physically 
more probable, and turbulence is indeed seen in the nu- 
merical calculations of Abramowicz, Igumenshchev & La- 
sota (1998). Recent work on an inhomogeneous solution to 
the accretion flow equations has shown that there is a sta- 
ble solution to the radiation pressure dominated SS disk 
if the majority of mass is in dense clumps which are opti- 
cally thick, while the remainder forms an optically thin(ish), 
X-ray hot plasma (Krolik 1998), with a physical mecha- 
nism for clumping provided by the photon bubble instability 
(Gammie 1998). However, one problem with these models is 
that they seem to predict that the covering fraction of cool 
clumps is close to unity for stellar mass black holes, in con- 
flict with the observations. 



9 CONCLUSIONS 

We detect significant relativistic smearing of the reflected 
spectral features in Cyg X-l, giving the first observational 
confirmation of an inner accretion disk in Cyg X-l. Our 
model calculates the total (line plus continuum) reflected 
spectrum for a given ionization state and iron abundance 
and then applies the relativistic smearing at a given incli- 
nation. These four parameters can be uniquely determined 
from X-ray spectra (assuming a given illuminating flux) 
with moderate energy resolution, and excellent signal-to- 
noise data extending beyond 9 keV. These conditions are 
satisfied in Cyg X-l for the EXOSAT GSPC and ASCA GIS 
spectra. The GINGA data are taken in a reduced spectral 
resolution mode, while the ASCA SIS data have insufficient 
signal-to-noise above 9 keV. 

For the ASCA GIS and EXOSAT GSPC data the re- 
flected spectrum is constrained to be from material inclined 
at > 45° which is not highly ionized. The low ionization state 
rules out Auger ionization as a mechanism for suppressing 
the line with respect to the reflected continuum as proposed 
by Ross et al (1996). The iron line is not weak with respect 
to the reflected continuum - we infer twice solar iron abun- 
dance - but the visibility of the line is dramatically reduced 
by relativistic smearing. 

The observed amount of relativistic smearing is not con- 
sistent with a disk extending all the way down to the last 
stable orbit around the black hole, but instead indicates that 
the disk truncates at some typical (but probably variable) 



© 0000 RAS, MNRAS 000, 000-000 



Relativistic Distortions in the X-ray Spectrum of Cyg X-l 13 



radius of ~ 10 — 30_R g , with a trend that larger transi- 
tion radii are associated with harder X-ray spectra. This 
mirrors the results from the hard spectra of the transient 
GBHC such as V404 Cyg and Nova Muscae (Zycki et al 
1997; 1998ab), suggesting that low state can be associated 
with a lack of optically thick disk at the innermost radii. 

The amplitude of the reflected component is signifi- 
cantly less than that expected from an isotropic source above 
a semi-infinite disk. Again this result can be understood in 
the context of an inner 'hole' in the optically thick material 
which is filled by the X-ray emitting plasma. Similar con- 
clusions about the geometry are also found from looking at 
the detailed continuum spectral shape of Cyg X-l, where 
the hardness of the X-ray spectrum implies that the hot X— 
ray plasma is spatially separate from most of the accretion 
disk, otherwise it intercepts too many soft photons from the 
disk and produces an X-ray spectrum which is significantly 
steeper than observed (Gierlihski et al 1997a; Dove et al 
1997; Poutanen et al 1997). 

The data then strongly indicate that low state GBHC 
accretion geometry is that of an inner X-ray hot, optically 
thin(ish), geometrically thick plasma, which cools by Comp- 
ton upscattering of soft photons from an external cool, op- 
tically thick, geometrically thin disk. Such a composite flow 
seems most likely to be a standard SS cool outer disk, with 
an inner ADAF or SLE flow. The SLE flow is probably 
thermally unstable, so for steady state solutions the ADAF 
seems to be preferred. However, we stress that the radiation 
pressure instability may put the disk in a regime where the 
flow is not in steady state, or is inhomogeneous. 
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